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ABSTRACT 



Context. Recent accurate abundance analyses of B-type main sequence stars in the solar vicinity has shown that abundances derived from 
these stellar objects are more homogeneous and metal-rich than previously thought. 

Aims. We investigate whether the inhomogeneity of abundances previously found in B-type stars in the OriOBl association is real (hence a 
signature of enrichment of the newly formed stars in an induced star formation scenario) or a consequence of intrinsic errors induced by the 
use of photometric indices to establish the stellar parameters prior to the abundance analysis. 

Methods. We obtained a new (improved) spectroscopic data set comprising 13 B-type stars in the various Ori OB 1 associations, and performed 
a detailed, self-consistent spectroscopic abundance analysis by means of the modern stellar atmosphere code fastwind. 

Results. We detect systematic errors in the stellar parameters determined previously which affect the derived abundances. Once these errors 
are accounted for, we find a high degree of homogeneity in the O and Si abundances for stars in the four Ori OB 1 subgroups. The derived 
abundances are in very good agreement with recent determinations in other B-type stars in the solar vicinity. We also compare our results with 
those obtained for the Sun during the epoch of the "solar crisis", and the Orion nebula. 
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;> 1. Introduction 



For many years, our knowledge about the chemical composi- 
tion of early-B main sequence stars in the solar vicinity has 
been characterized by two main results: (i) the derived abun- 
dances seemed to be highly inhomogeneous (with a dispersion 
of up to 0.5 dex), and (ii) the mean values indicated lower abun- 
dances than the standard ( Grevesse & Sauval , 1998I) set of solar 
abund ances (see reviews by He rrero & Lennoni l2004t iMorel , 
2009). These results were not very encouraging, because the 



younger than the Sun, nearby OB-type stars are expected to be 
slightly metal-rich). 

Some recent results have began to change this sit- 
uation. The solar oxygen abundance t raditionally consid- 
ered as a cosmic abundan ce reference (iGrevesse & SauvalL 
1998) was reviewed by lAsplund et al.1 J2004 ~ 



who ob- 

0.17 dex lower than the stan- 
beginning of what has been 
crisis": between 2004 and 



inhomogeneity of stellar abundances contradicted with the ho- 
mogeneity in oxygen abundanc e found from studying the local 
diffuse interstel lar medium (e.g. lMever. Jura. & Cardelilll998l: 



Cartledge et al. , 2006). On the other hand, chemical evolution 



tained log (0/H)=8.66 dex, 
dard value. This was the 
called the epoch of the"solar 

2008 , several studies by di f ferent authors dAvres et al 



20061: ISocas-Navarro & Nortonl l2007t lAllende Prietol 



Caffau et al 



20081: lAvresl 1 20081 TMelendez & Asplund . 



2008; 



2008 



Centeno & Socas-Navar ro, 2008) presented solar oxygen abun- 



models of the Ga l axy (e .g. IChiappini. Romano. & Matteucci , 
l2003t l ICarigi et al.ll2005l) predict a small enrichment of the ISM 
in metals during the lifetime of the Sun (i.e. because they are 



dances derived by means of differe nt approaches. The calcu- 
lated values range from 8.63 dex (fSocas-Navarro & Nortonl 



20071) to 8.86 dex (ICenteno & Socas-NavarroL 120081) . The de- 
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* Based on observations made with the Nordic Optical Telescope, 
operated jointly on the island of La Palma by Denmark, Finland, 
Iceland, Norway, and Sweden, in the Spanish Observatorio del Roque 
de los Muchachos of the Instituto de Astrofisica de Canarias. 



bate abou t its actual value is st i ll ope n. 

IPrzvbilla. Nieva. & Butlet ( 2008 ) have recently analyzed a 
representative sample of six unevolved early B-type stars in 
nearby OB associations and the field, and found a very nar- 
row distribution of abundances, with mean values that are more 
metal-rich compared to previous works (e.g. log (0/H)=8.76 
dex, a value that is within the range of solar abundances cal- 
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culated during the "solar crisis")- These authors indicate the 
importance of properly determining the atmospheric parame- 
ters and using robust model atoms to avoi d systematic errors 
in the abundance determination. (See also INieva & Przv billa. 
2009, for a summary of the main sources of systematic errors 
aff ecting the abundance analy ses of B-type stars.) The study 
by IPrzybilla. Nieva. & Butlerl show that the chemical inhomo- 
geneity previously found for B-type stars in the solar vicinity 
may be spurious and an artificial effect of those systematic er- 
rors. It also reinforces the importance of self-consistent spec- 
troscopic abundance analyses (i.e., the stellar parameters and 
the metal abundances are determined exclusively from spectro- 
scopic diagnostics by using the same set of stellar atmosphere 
models). Less accurate photometric T e g estimates must be han- 
dle with care (or better avoided wheneve r possible!) in the 
abundance analysis of B-type stars (see also lNieva & Przvbilla , 



2008). 



The Orion complex, containing the Orion molecular cloud 
and the Orion OBI (Ori OBI) association, is one of the most 
massive active star-form ing regions in the 1 kpc centered on 
the Sun. iBlaauwl (|l964) divided Ori OBI into four subgroups 
of stars — namely la , lb, Ic, and Id — h aving different locations 



in the sky and ages. Brown et al.1 (119941) derived mean ages of 



11.4+1.9, 1.7+1.10,4.6+2, and <1 Myrfor subgroups la to Id 
respectively. The youngest subgroup Ori OB 1 Id is associated 
with the Orion nebula (M 42), the most studied H n region and 
the closest ionized nebula to the Sun in which a high accuracy 
abundance analysis can be performed. 

The correlation between the ages of the stellar subgroups, 
their location, and the large scale structures in the interstel- 
lar medium around Orion OBI have been interpreted as fea- 
tures of sequential star formation and Type- I I supernovae 



dRevnolds & Qgdenl fl979) ICowie et all 1 19791 iBrown et al. 



1994 . ICunha & Lambertl (U9921 Il994h obtained C, N, O, Si 



and Fe abundances of 18 B-type main sequence stars from 
the four subgroups comprising the Ori OBI association. They 
found a range in oxygen abundances^ of ~0.4 dex, with the 
highest values corresponding to the stars in the youngest (Id 
and some Ic) subgroups. In this case, the inhomogeneity in 
stellar abundances (mainly oxygen and silicon) seemed to be 
real and coherent with a scenario of induced star formation in 
which the new generation of stars are formed from interstellar 
material contamina ted by Type-II supe rnovae ejecta. 

The study by ICunha & Lambert! was based on a pho- 
tometric estimatio n of T f ff and the fitting of the Hy line 
computed from Kurucz ( 1979t) lte model atmospheres to 
the observed one to de rive logg. In a more recent work, 
Simon-Diaz et al. d2006l) used one of the new generation 



of nlte, line blanke t ed, model atmosp h ere co des (fastwind, 



Santolava-Revetal., 1997; Pulsetal., 2005) and a self 



1 iBriceno etall d2005h determined an age for Ori OB lb ~4 - 6 Myr 
by studying its low mass, young stellar population; this value is more 
consistent with what would be expected from the presence of the 
evolved blue supergiant e Ori in this subg ro up. 

2 Other authors dGies & Lambertl Il992l ; iKilianl 1 19921 ; 
iGummersbach et all Il998h obtained a similar range in oxygen 
abundances from the analysis of smaller samples of B-type stars in 
Orion OBI. 



consistent spectroscopic approach to derive the stellar param- 
eters and oxygen abundances for the three B0. 5 V stars in 
Ori OB Id (the youngest and, supposedly, more metal-rich sub- 
group). The resulting stellar parameters were somewhat differ- 
ent and, more important, the derived abund ances were system- 
aticall y lower than the previous values by ICunha & Lambertl 
d 19941) by ~ 0.2-0.3 dex. 

This result motivated us to review the chemical composi- 
tion of the other B-type stars in Ori OB 1 to investigate whether 
the inhomogeneity of abundances previously found is real or 
a consequence of intrinsic errors induced by the use of pho- 
tometric indices to establish the stellar parameters prior to 
the abundance analysis. To this aim, we obtained a totally 
new (improved) observational data set and performed a self- 
consistent abundance analysis o f 1 3 of the stars considered by 
Cunha & Lambertl d 19921 1 1994 . 

We used fastwind to derive the stellar parameters, oxygen, 
and silicon abundances. The observational data set is described 
in Sect.|2] The whole spectroscopic analysis is presented in 
Sect.|3] Then, we compare our stellar parameters and abun- 
dance with results from previous works (Sect.[4]>. The homo- 
geneity of stellar abundances in Ori OBI and its comparison 
with other B-type stars in the solar vicinity determinations, the 
Sun, and the Orion nebula is discussed in Sect. [5] The main 
conclusions of this work are summarized in Sect.|6] 



2. Observational data set 

The observations used here were carried out with the FIES 
cross-dispersed, high-resolution echelle spectrograph attached 
to the NOT2.5m telescope at El Roque de los Muchachos 
observatory on La Palma (Islas Canarias, Spain) on 5-8 
November 2008. The medium-resolution mode (R=46000, 
5/1=0.03 A/pix) was selected, and the entire spectral range 
3700-7300 A was covered without gaps in a single fixed set- 
ting. A sample of 14 stellar candidates was observed, after be- 
ing selected from the list of stars in the Or i OBI association 
analyzed by lCunha & Lambertl d 199211 1994 . These are early- 
B type, non-evolved, stellar objects (B0 V-B2 V) with low pro- 
jected rotational velocities (vsin/<60 kms ). The signal-to- 
noise ratio achieved for all the spectra was always above 250. 
The list of observed stars is presented in TableQ] Apart from the 
new set of stars, we re-observed HD 37020 and HP 37042, tw o 
of the three B-type stars^l analyzed in lSimon-Dfaz et aT (2006). 

The spectra were reduced with the FIEStoofl software in 
advanced mode. The FIEStool pipeline provided wavelength 
calibrated, blaze-corrected, order-merged spectra of high qual- 
ity. These spectra were then normalized with our own devel- 
oped idl programs. An example of the fies@not spectra for 
three of the observed stars is presented in FigQ] where the Si n- 
iv and O n lines used for the abundance analysis are also indi- 
cated. 



3 HD 37023 was also re-observed, but it was not included in the 
analysis because we found in the new spectrum that the star is is actu- 
ally a binary system (SB2), hence not optimal for an abundance anal- 
ysis. 

4 http://www.not.iac.es/instruments/fies/fiestool/FIEStool.html 
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Fig. 1 Example of fies@not spectra for three of the observed stars. The complete atlas is available in the electronic version of 
the paper (Fig.fTOil. The Sin-iv and On lines used for the abundance analysis are indicated as solid and dashed vertical lines, 
respectively. 



3. Spectroscopic analysis 

The analyses were performed following a self-consistent 
spectroscopic approach with the spherically extended, 
nlte, line-blanketed stellar atmosphere code fastwind 



(Santolav a-Rev et all 11997; Puis et al.L 120051) . Basically, the 



stellar parameters were derived by comparing the observed H 
Balmer line profiles and the ratio of Si m-iv and/or Si ii-m line 
equivalent widths with the output from a grid of fastwind mod- 
els. Whenever possible, the Hei-n ionizing equilibrium was 
also considered. Then, the same grid of models was used to 
derive the stellar abundances by means of the curve-of-growth 
method. 

3.1. Projected rotational velocities 

The projected rotational velocities (vsin Q were obtained 



by applying the Fou r ier method (Gray] Il976t see also 



Simon-Diaz & Herrerol I2007L for its application to OB-type 



stars) to the Sim/14552 line. The 0.03 A/pix resolution of the 
FIES spectra implies that the lowest v sin i that could be de- 
tected is 2-3 kms~'; however, for those cases with vsin/< 10- 
15 km s _1 , identifying of the first zero of the Fourier transform 
was difficult due to the effect of the noise (an maybe the micro- 
turbulence). 

In many cases, an additional extra Gaussian-type broad- 
ening (& g ) was needed to properly fit the line profile for the 
derived value of v sin i. This extra-broadening account for the 
microturbulence (and maybe the macroturbulence), also affect- 



Table 1 List of B-type stars from Ori OBI considered in this 
study, ordered by spectral type. 

HD HR Name SpT V Subgroup 



HD 36512 


HR1855 


y Ori 


B0V 


4.62 


Ic 


HD 37020 


HR1893 


9 [ Ori A 


B0.5 V 


6.71 


Id 


HD 36960 


HR1887 




B0.5V 


4.78 


Ic 


HD 37042 




<9 2 0riB 


B0.7V 


6.02 


Id 


HD 36591 


HR1861 




Bl V 


5.34 


lb 


HD 36959 


HR1886 




Bl V 


5.67 


Ic 


HD 35299 


HR1781 




B1.5 V 


5.69 


la 


HD 37744 


HR 1950 




B1.5 V 


6.22 


lb 


HD 36285 


HR1840 




B2V 


6.33 


Ic 


HD 36629 






B2V 


7.65 


Ic 


HD 35039 


HR 1765 




B2V 


4.72 


la 


HD 36430 


HR1848 




B2V 


6.23 


Ic 


HD 35912 


HR1820 




B2V 


6.41 


la 



ing the shape of the line. The corresponding derived values are 
summarized in Table [2] 



3.2. Measurement of equivalent widths and 
identification of problematic lines 

The strategy we followed in our analyses is based on the equiv- 
alent widths (EW) of metal lines. Therefore, proper identifi- 
cation of the lines of interest, along with an accurate measure- 
ment of their EWs (also including the associated uncertainties), 



4 



S. Simon-Diaz: O and Si abundances in B-type stars in Ori OB 1 



Table 2 Projected rotational velocities derived for the studied 
stars. 



HD 


v sin i" 


e; 


HD 


v sin i 


®i 


HD36512 


15 


30 


HD 37744 


37 


00 


HD 37020 


55 


00 


HD 36285 


<10 


20 


HD 36960 


23 


40 


HD 36629 


<10 


15 


HD 37042 


31 


00 


HD 35039 


<12 


15 


HD 36591 


08 


19 


HD 36430 


<15 


25 


HD 36959 


12 


12 


HD 35912 


<12 


20 


HD 35299 


<08 


12 









"Both quantities in km s 



is a very important step. We have developed IDL routines auto- 
matically identifying metal lines in the spectra, measuring the 
EWs and their uncertainties, and detecting the possibility that 
other lines affect these measurements. To this aim we compiled 
a list of C, N, O, Si, Mg, S, Ne, and Ar lines, extracted from 
the atomicjfl line list v2.05. For a given line, the program per- 
forms a multi-Gaussian fit of the observed line profile account- 
ing for all the lines expected to be present in a certain spectral 
range (/lo±2max[vsin//lo/c, 0.5/to/R]) around the wavelength 
indicated in the line list. The uncertainty in the EW measure- 
ment is obtained by assuming the location of the local contin- 
uum at +1/SNR and, in those cases in which the line is isolated, 
comparing the value obtained by means of the Gaussian fitting 
with the value derived by integrating the line. 

The high quality of our spectra allowed us to identify and 
measure the EWs of up to 27 Sin-iv lines and 47 On lines. 
Some of the lines were labeled as problematic because of the 
presence of one or more lines from other elements (e.g. the 
O n/14673.73 line may be affected by the C ni/14673.95 in same 
cases; similarly occurs for the Oii/14641.81 line, coincident 
with the N m /14641 .85 line). These lines were treated with spe- 
cial care in the abundance analysis, since they may be giving 
wrong values for the abundance. 

Generally, the Gaussian fit provides reliable results for the 
EWs; however, in those cases in which the vsin/ of the star 
is above 30-35 kms -1 , the use of a Gaussian to fit the profile 
may result in an under or overestimation of the EW (in a few 
percent), depending on the line strength. For those cases, the 
EW resulting from the integration of the observed line profile 
was preferred. 

3.3. A grid of HHeOSi fastwind models 

For this study, we constructed a grid of fastwind models with 
Teft and log g ranging from 17 to 36 kK (1 kK steps) and 3.7 to 
4.3 dex (0.1 dex steps). As the studied stars are not expected 
to be evolved, the He abundance was fixed to 0.09 dex. In 
addition, since fastwind is a spherically extended code, the 
radius and other wind parameters need to be indicated, and 
are grouped in the Q-parameter. We fixed this parameter to 



logQ=-15 as a representative value for which the wind ef- 
fect over the optical spectrum is practically negligible. (Ha and 
He 11/14686 show no sign of wind contamination.) The metallic- 
ity was assumed to be sola r (following the set of abundances by 
iGrevesse & Sauvallll998l) . 

For each pair of stellar parameters, a sub grid of models 
varying the microturbulence (ft=l> 3, 5, 7, 9 kms^ 1 ), the Si 
abundance (e si =-5.10, -4.80, -4.50, -4.20 dex), and the O 
abundance (eo=^-00, -3.65, -3.30, -2.95 dex) was calcu- 
lated. The O and Si atomic models used for the grid came 
mainly from lBecker & Butler ( 1988 , 1990)- However, two up- 
dates were consider ed: (a) a extended Sin model atom (see 
iTrundle et all l2004h . and (b) the most recent log gf values in- 
dicated in the atomic line list v2.05 for the formal solution cal- 
culations. 

The final grid consists of 20x7x4 (=560) models and a total 
of 2800 formal solutions (5 microturbulence values per model). 
It includes line profiles and EWs for H , He i-n, Si n-iv and O n 
lines, along with the spectral energy distribution for each set of 
stellar parameters. 

3.4. Determining stellar parameters 



The use of the Sim-iv and/or Sin-m ionization equilib- 
rium, along with the H Balmer lines, for determining the 
stellar parameters of early B-type stars is a longs t anding 
method described e l sewhere (see e.g. Kilian et all 1991 



http://www.pa.uky.edu/~peter/newpage/ 



l Urbaneia etall 120051: ICrowther et all 120061: iMarkova & Puis . 
2008]). The ratios W(SiiV/l4116)/£W(Siin/l4552) and/or 
W(Siii/l4128)/£W(Siiii/l4552), depending on the tempera- 
ture of the star have been traditionally used as T e jf indicators. 
This decision has been probably motivated by the fact that these 
are the strongest, unblended lines in the spectral range com- 
monly observed for the stellar abundance analysis (i.e. ~ 4000- 
5000 A). 

We derived the stellar parameters, along with the Si abun- 
dance, for all stars in our sample using this methodology. We 
initially considered that these Si line ratios determine the ef- 
fective temperatures (columns 3 and 4 in Table [3). However, 
motivated by a couple of problems found when deriving the 
stellar parameters in the cooler objects (see discussion below), 
we decide to include the £W(Sin/l6347)/£ , W(Siiii/14552) ra- 
tio as a temperature indicator (column 5). The measured values 
are indicated in brackets in the corresponding columns. As ex- 
pected, the Siiv/Sim ratio decreases when we move to later 
spectral types, and the Si n/Si m behaves in the opposite way. 
For three of the stars (HD 36959, HD 35299, and HD 37744), 
lines from the three ions are clearly and simultaneously present 
in the spectra. 

Four parameters need to be determined at the same time in 
an iterative way: T e ff, logg, £ t (Si), and es;. First, we use the 
Si line ratios indicated above and the wings of the H Balmer 
lines (fixing £(Si), and e^) to obtain an initial guess for T e ff and 
log g. Then we apply the curve of growth method to a proper set 
of Si n-rv lines to iteratively obtain final values for the four pa- 
rameters (a detailed description of the used lines and the results 
of the Si abundance analysis is presented below). Normally, the 
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Table 3 Final results of the HSi analysis: stellar parameters and Si abundances. Detailed results (line-by-line) of the Si abundance 
analysis are presented in Tables [6lto[T8lin the electronic version of the paper. All T e ff values correspond to the \ogg indicated in 
the 8th column. 







(EW4/EW,)" 


(EW 2A /EW 3 ) 


(EW 2 b/EW]) 




Adopted 








TW+500K 




SpT 


T „ TV— TFT 

± gff — ijl 1 V 111 


T rr TT-TTT 

1. gff — ui 11 ill 


T ,,. Ql TT _ TTT 

± gff kjl 11 111 


Tgff — He i~ ii 


T e ff 


log,? 


A ("Si") 


ft.- -+- 


Aft-Trr £S\ 
^tSiL » St J 


^tSi V^ 1 


HD 36512 


B0 V 


(1.42±0.06) 
33700±200 






34000+500 


33700 


4.2 


4.3+0.7 


7.49+ 


[0.07, 0.05] 


±0.12 


HD 37020 


B0.5 V 


(0.69±0.09) 
30500±600 






30000+500 


30500 


4.2 


0.5+0.5 


7.47+ 


[0.10, 0.04] 


±0.06 


HD 36960 


B0.5 V 


(0.60±0.04) 
28900±300 






29000+500 


28900 


3.9 


5.4+0.6 


7.53+ 


[0.02, 0.06] 


±0.06 


HD 37042 


B0.7 V 


(0.64±0.07) 
29700±400 






29500+500 


29700 


4.2 


1.4+0.3 


7.55+ 


[0.03, 0.04] 


±0.05 


HD 36591 


Bl V 


(0.39±0.01) 
27200±100 








27200 


4.1 


1.3+0.3 


7.53+ 


[0.06, 0.03] 


±0.01 


HD 36959 


Bl V 


(0.25±0.02) 
25900±300 


(0.12*0.02) 
25900+100 






25900 


4.2 


0.0+1.0 


7.50+ 


[0.05, 0.07] 


+0.02 


HD 37744 


B1.5 V 


(0.16±0.03) 
23900±600 


(0.12±0.01) 
25700+100 


(0.12±0.04) 

23600+600 




23800 


4.1 


0.5+0.5 


7.54+ 


[0.06, 0.04] 


+0.05 


HD 35299 


B1.5V 


(0.14±0.02) 

23900±300 


(0.20±0.01) 
24900±300 


(0.17*0.03) 
23700±300 




23700 


4.2 


0.5+0.5 


7.50+ 


[0.08, 0.02] 


+0.06 


HD 36285 


B2V 




(0.26±0.01) 

23900±200 


(0.41 ±0.03) 

20600±200 




20600 


4.0 


1.7+0.5 


7.49+ 


[0.06, 0.05] 


+0.11 


HD 35039 


B2V 




(0.36±0.02) 
22200±200 


(0.44*0.04) 

19900±200 




19800 


3.7 


3.3+1.0 


7.52+ 


[0.06, 0.08] 


+0.14 


HD 36629 


B2V 




(0.39±0.02) 

22800±300 


(0.53*0.02) 

20000±100 




20000 


4.1 


1.0+0.5 


7.54+ 


[0.04, 0.05] 


+0.13 


HD 36430 


B2V 




(0.65±0.06) 

21000±300 


(0.96*0.09) 

18600±200 




18600 


4.1 


3.5+1.0 


7.47+ 


[0.08, 0.07] 


+0.13 


HD 35912 


B2V 




(0.73±0.05) 

20400+300 


(0.92*0.05) 

18500±150 




18500 


4.0 


3.2+0.5 


7.48+ 


[0.07, 0.04] 


+0.13 



a EW 4 : Siiv4116; EW 3 : Sim 4552; EW 2A : Sin4130; EW 2B : Sin 6371 



final values of T e ft and logg are quite close to the initial values, 
since the mentioned line ratios are only slightly dependent on 
e si and 



3.4.1 . Hotter objects (T eff > 27000 K) 

In the hotter objects, Si n lines are not present in the spectra. 
The EW(Si iv A41 16)/£W(Si m ,14552) ratio is thus used to ob- 
tain the initial guess values of T e ff for each of the log g val- 
ues considered in the grid. Then, the full set of reliable Si m-iv 
lines, along with the H Balmer line profiles, is used for fine 
determination of the four parameters indicated above (see an 
example in the upper panel of Fig. [2] and in Fig. [3j. The un- 
certainties in the EW measurements of the Si lines have been 
taken into account for establishing the uncertainties associated 
with the derived stellar parameters. Generally, the T e ff and log g 
can be determined with an accuracy better than 500-600 K and 
0. 1 dex, respectively. The final results of the analysis are pre- 
sented in Table[3] 

For main sequence stars with T e fj > 28000 K the He i-n ion 



ization equilibrium can also b e used (see e.g. iHerrero et al 



19921: [Simon-Dfazetail [2006). (Below this temperature He n 



lines are too faint or not present in the spectrum.) We could de- 
termine the stellar parameters in this way for four of the stars 
in our sample. A n example of this type of analysis has already 
been presented in lSimon-Dfaz et al. (2006), so we only present 
here the corresponding results (see column 6 in Table [3]). In 
general, there is good agreement between the T e tj determined 
through the Si m-iv and the He i-ii ionization balance (with dif- 
ferences in T e ff not larger than ~500 K. 



3.4.2. Cooler objects 

For the cooler objects, the Si n/Si m must be used. We initially 
considered the £'W(Sin^4130)/£ , W(SiiiU4552) ratio for ob- 
taining a first estimation of the stellar parameters; however, we 
found two facts that warned us of a possible problem with the 
Siii/14130 (or the Siii/14128) line. First, it was not possible 
to properly fit the H Balmer lines for any of the (T e ff, logg)- 
pairs indicated by the Si n /14130/Si m /14552 line ratio, because 
the core of the H synthetic lines were somewhat narrower than 
the observed ones (even if high values of log g are considered). 
For these temperatures, the cores of the H Balmer lines be- 
gin to be sensitive to changes in T e ff. For all the cases stud- 
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ied, these lines require somewhat lower effective temperatures 
to properly fit their cores. Second, the egj-EW diagnostic dia- 
grams show two sets of Si n lines giving different abundances 
by ~0.2 dex (see lower panel in Fig.|2j». Curiously, each subset 
of lines corresponds to transitions with very different energy 
levels (Si ii/t/14128, 4130, 5041, 5056, in one hand, coming 
from higher energy levels; Siii/L16347, 6371, 3856, 3862, on 
the other, comming from lower energy levels; we suggest the 
reader have a Grotrian diagram in hand). This different behav- 
ior of the various lines warned us of possible problems with the 
Sin atomic model and of using Siii/14128 line alone to estab- 
lish the stellar temperature. 

Columns 4 and 5 in Table [3] show the different ef- 
fective temperatures obtained depending on the Sin line 
that is used (Siii/14128, Siii/16371). Temperatures given by 
the W(Siii/l6371)/£W(Siiii/l4552) ratio are systematically 
lower by up to ~2000- 3000 K. In addition, a proper fit of the 
cores of the H Balmer lines can be achieved with these new 
values of the temperature. 

3.4.3. HD 36595, HD 37744, and HD 35299 

These three stars show Sin-m-iv lines in their spectra, so 
can contribute with new decisive clues to the problem men- 
tioned above. In this case we count with three T e ff indi- 
cators. As can be noticed from inspection of Table [3] the 
T eff indicated by the £W(SiiU6371)/£W(Sim/l4552) ratio 
results in better agreement than the other Sin/Sim ratio for 
two of the three cases. For the hotter object (HD 36959), the 
Siii/16371 is too faint to be measured. Although the T e ff in- 
dicated by the EW(Sin ,14 130)/W(Sim ,14552) ratio agrees 
with the Si iv/Si m diagnostic in this object, this is not the case 
for the other objects. 

In Fig. [2] (middle panel) we show the esi-EW diagnostic di- 
agram for HD 37744. The stellar parameters considered in this 
plot are those indicated by the Si iv/Si m ratio. The Si n AA637 1 , 
3856 lines fit the other Sim-iv lines perfectly, but not the 
Si H/l/14128,4130 lines. 



3.4.4. Concluding 

We found some indications of a problem in the Si n model atom 
that led to bad modeling of the Sin/U4128, 4130, 5041, 5056 
lines. Fortunately, our observations also include other Si n lines 
that seem to behave properly. A review of the Si n model atom 
is needed (and could be tested with the type of detailed analysis 
presented here), but in the meantime, several arguments allow 
us to trusjl the Si H/U6347, 6371, 3856, 3862 set of lines for 
the stellar parameter determination and for abundance analysis 
of the cooler objects in our sample: (a) the coherence of results 
between Si iv-m and Si m-ii in terms of T e ff ; (b) the good fit of 
the H Balmer lines for the cooler stars when the lower T e ff is 
considered; (c) the coherence of results in Si and O abundances 



HD 36960 - T Bff =28900 K, logg = 3.9, {-, = 5.4 



6 This hypothesis is also supported by a comparison of T e g obtained 
for the cooler objects in our sample from other spectroscopic diagnos- 
tic based on O ii/i and C m/n ratios of lines (F. Nieva, private commu- 
nication). 
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Fig. 2 Silicon abundance vs. EW diagnostic diagrams for three 
of the analyzed stars (with representative effective tempera- 
tures). 



that is obtained for all the stars in our sample when this solution 
is adopted (see below). 



3.4.5. Uncertainties 

Columns 3 and 5 indicate the uncertainties in T e ff obtained 
from considering errors in the EWs measurement of Sin-iv 
lines. These uncertainties, obtained by assuming a fixed gravity 
are ~ 1 00-600 K. Gravity can be normally established with an 
accuracy better 0. 1 dex. Uncertainties from both quantities are 
correlated, e.g. a positive variation of log g of 0.1 dex needs to 
be compensated by an increase of T e ff ~1 00-300 K to guaran- 
tee again the Si ionization balance. Although in many cases the 
formal errors in temperature obtained from the propagation of 
EWs errors are smaller than 500 K, our experience warns us to 
be conservative (since other parameters can also slightly affect 
the derived temperatures, such as the considered microturbu- 
lence). Therefore, we adopt 500 K and 0.1 dex as characteristic 
uncertainties in T e ff and log g, respectively, from our analyses. 
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Fig. 3 Example of fitting the H Balmer lines to determine logg. 
Two fastwind models are compared with the observed spec- 
trum (red line) of HD 36960: in black, model with the stellar 
parameters resulting from the fastwind HSi analysis (28900, 
3.9); in blue, model w ith the stellar parameters derived by 
Cunha & Lambertl(ll992h . for comparison (28900, 4.3). 



3.5. Silicon and oxygen stellar abundances 

We applied the curve-of-growth method to derive Si and O 
abundances. This method considers a grid of models for a given 
set of stellar parameters (T e jf and log g) in which the microtur- 
bulence and the abundance of the element to be studied are var- 
ied. An abundance is obtained for the various values of micro- 
turbulence for each of the considered lines (given the measured 
EW of the lines in the observed spectrum). Then, the final abun- 
dance is given by that microturbulence that results in all lines 
giving the sam e abundance . More details on this meth od can be 
found in e.g. lKilianl (Il992h or lSimon-Dfaz et al.1 d2006l) . 

Our preference for this methodology in performing the 
abundance analysis is that we find the curve-of-growth method 
very powerful for identifying problematic lines (as shown e.g. 
in the previous section), which can affect the final abundance 
determination, and provide precise estimations of the uncer- 
tainties associated with the dispersion of line-to-line abun- 
dances, microturbulence, and the stellar parameters. 

3.5.1. Silicon abundances 

As mentioned, silicon abundance is obtained at the 
same time as to the stellar parameters. Once a first 
estimation of the stellar parameters was obtained 
through the £W(Sirvi4116)/£W(Siiii/l4552) and/or 
£W(Sin ,1637 l)/£W(Sim ,14552) line ratios, the curve- 
of-growth method was applied to a proper set of Si lines to 
derive the Si abundance, together with the final values of the 
stellar parameters and a microturbulence. 

Initially, we included all the measured lines in the analy- 
sis, but identified some problematic lines. By problematic lines 
we mean lines giving abundances that are too high or too low 
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Fig. 4 Silicon abundance vs. EW diagnostic diagrams for 
HD 36960. All the observed Si ni-iv lines are included in the 
plot. Uncertainties in the individual line abundances (propa- 
gated from the errors in the measured EWs) are indicated as 
vertical lines 



compared with the mean value provide by a set of lines initially 
considered as reliable. These problematic lines were removed 
in the final analysis. To illustrate the procedure, we consider 
the case of HD 36960. The esi-EW diagram with whole set of 
measured Si lines (Si m and Si iv lines, in this case) is shown 
in Fig. |U As the initial set of reliable lines we consider the 
Sim/U4552, 4567, 4574 triplet and the Si iv ,141 16 line. The 
other Si iv lines provide similar abundances to the first set of 
reliable lines, except Si iv ,16701. The situation is less encour- 
aging for the other Si m lines, since all of them lie below the 
mean abundance value from reliable lines. We thus label these 
lines as problematic. The same procedure was followed for the 
selection of reliable Sin lines (see Sect. 13.41 ). 

A similar behavior is found for the problematic lines in 
all the analyzed stars. Since the EW of these lines are accu- 
rately measured and no lines from other elements are expected 
to be affecting therrQ, we argue that the discrepancy may be 
related to the definition of the model atom. For example, the 
Sim/U4813, 4819, 4829 triplet is known to give different re- 
sults than the Si m AA4552, 4567, 4574 tripl et because of the 
boun dary problems of the Si m model atom dBecker & Butlei , 
1990). Normally, there is an explanation for the bad behavior of 
the problematic lines; therefore, the exclusion of lines from the 
analysis is not arbitrary (e.g. lines from the same multiplets be- 
have in the same way and are normally excluded all together). 

It is important to notice that the analysis with the whole 
set of lines (including the problematic ones) could lead to an 
incorrect determination of the microturbulence, hence of the 
final value of Si abundance. 



7 This is not the case for Si iv ,1,14089, 463 1 , and Si n ,1,16347, 3856; 
these lines may be blended with On ,14089, Nn ,14631, Mgn ,16347, 
and O n/13856, respectively, in some cases. The high resolution of the 
fies@not spectra normally allows separating both line contributions; 
however, results with these lines are always treated with care. 
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Final results of the Si abundance analysi^] are indicated in 
Table [3] (columns 9 and 10). Three sources of errors must be 
considered for estimating of the final uncertainty: (1) the dis- 
persion in the line-by-line abundances, (2) the uncertainty as- 
sociated with the microturbulence, (indicated in Table [3] and 
Tables l6lfT8l only available in the electronic version), and (3) 
the contribution of uncertainties in the stellar parameters. These 
can be added quadratically to obtain the total uncertainty. 

For illustrative purposes in the last column of Table [3] we 
indicate the effect of a change of +500 K in T e jf on the Si 
abundance if this was derived exclusively by using Si in lines^l 
The minimum effect occurs for T e ff~ 27000 K and increases 
towards higher and lower temperatures (note the change in the 
sign of the uncertainties), reaching values up to 0.12 dex. This 
behavior is a consequence of the dependence of the EW of the 
Sim lines with temperature. The maximum EW is achieved 
around 27000 K (so the abundance is quite insensitive to small 
changes in T e ff ), and decreases towards lower and higher tem- 
peratures (also increasing its sensitivity to T e ff variations). 
When lines from two different ions are used, the final Si abun- 
dance can be constrained with better accuracy, since the EW of 
lines from different ions behaves in a opposite way for a given 
T e jf. As a consequence in this case, the uncertainty associated 
with the stellar parameters is always negligible compared with 
the dispersion in the line-by-line abundances. 

3.5.2. Oxygen abundances 

Up to 47 On lines were identified in the observed spectral 
range. However, not all the lines were finally used for the O 
abundance determination. The selection of the final set of lines 
used for the analysis was based on a detailed analysis by multi- 
plets (following similar criteria to the case of Si, but in this case 
we studied the behavior of lines resulting from different multi- 
pl ets). We found a similar b ehavior of lines to the one described 
in lSimon-Dfaz et all d2006h . 

The derived O abundances are indicated in Table|H and the 
results of the line-by-line analyses in Tables I61TT81 (in the elec- 
tronic version). For this element, only lines from one ionization 
state were available, so we could not test whether the corre- 
sponding ionization equilibrium (Om/ii or Oii/i) is achieved 
for the considered stellar parameters. We adopt for the oxy- 
gen abundance analysis the same stellar parameters as in the Si 
analysis. 

We found (as in many previous works) that the microturbu- 
lence derived from the O n lines (^t(O)) differs from the Si anal- 
ysis (£ t (Si)), and a somewhat larger microturbulence is derived. 
Some authors assume the £ t (Si) value (or a mean value of the 
microturbulences obtained for the various elements analyzed) 
to perform the oxygen abundance analysis. In our opinion, this 
can lead to significant systematic errors in the analysis. Since 
this is an ad-hoc parameter that is still not well understood (see, 



however, ICantiello et al. , 2009), we adopt the microturbulence 
derived from the oxygen analysis itself, for consistency. 

In fact, determination of the microturbulence value that 
will be adopted in the final steps of the abundance determi- 
nation is an important task. Unidentified problems in the O n 
line modeling or bad measurements of the corresponding EWs 
(due to blends, noise, or a bad placement of the continuum) 
can enormously affect the £ t value that produces a zero slope 
in the ep-EW diagr ams. A detailed analysis by multiplets (see 
Simon- Dial [2005) can help identify problematic lines and to 
better decide on the final microturbulence to be adopted. 

Table H also indicates the uncertainties associated with er- 
rors in the line-to-line abundance dispersion, the microtur- 
bulence, and the stellar parameters. In addition, last column 
shows the derived abundances if the effective temperature is 
varied +500 K. (Although the exact values are not shown here, 
the contribution of the log g uncertainty to the oxygen abun- 
dance can be considered negligible in comparison with the T e ff 
contribution.) The derived oxygen abundance is very sensitive 
to changes in T e g for the cooler and hotter objects, and there 
is a change in the behavior of the oxygen abundance with T e ff 
around 27000 K. This behavior is similar to the one illustrated 
in Table[3] 

4. Comparison with previous works 

Two of th e stars included in t his an alysis have already been an- 
alyzed in ISimon-Dfaz et al. | d2006l) . These are HD 37020 and 
HD 37042. We observed the two stars again to have spectra 
with the same characteristics as the other stars in the sam- 
ple. The stellar parameters presented in Table [2] correspond to 
the analysis of the new spectra. The new analysis resulted in 
slightly higher effective temperatures (but within the errors) 
and gravities -0.1-0.2 dex larger. Several factors produced 
this difference in the derived stellar parameters. First, we based 
our T e ff - log g determination on the HSi analysis, instead of us- 
ing the H and He i-ii lines. The best T e ff - log g pair reproducing 
the Si iii-iv ionization equilibrium and the wings of the H lines 
simultaneously is the one indicated in Table[3] A lower log g re- 
quested an effective temperature somewhat lower (which was 
not fitting the Hen lines). In addition, the better quality of the 
new spectra allowed us to better constrain the gravity of the 
stars. The wings of the H lines are not very sensitive to changes 
of ~ 0.1-0.2 dex in logg in this range of stellar parameters. In 
Simon-Diaz et al.. we based our decision about the best solu- 



tion on the faint Hen/14541 line; however, it is better to rely 
on the HSi criterion, which is more sensitive in this range of 
stellar parameters. In fact, we find that the new T e ft, - log g pair 
also fits nicely the He i-ii lines. We thus prefer this last solution. 

We also obtain ~0. 05-0. 10 dex higher oxygen abundances. 
Part of this difference is caused by the change in stellar 
parameters; in addition, we measured sligthly larger EWs for 
the On lines in both stars. Note, however, that the old and new 
abundances agree when taking the error bars into account. The 



8 Results of the line-by-line analysis are presented in Tables 1 61181 
(only available in the electronic version). 

9 A similar test assuming a change of 0. 1 dex in log g indicates neg- 
ligible effects when compared with the T c ff contribution. 



discrep ancy for these two objects between ISimon-Dfaz et al 



(2006) and these new studies serves as an example of the 
effect of intrinsic uncertainties in determination of stellar 
abundances. Following a similar methodology the same author 
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Table 4 Final results of the O abundance analysis. Detailed results (line-by-line) of the O abundance analysis are presented in 
Tables|6lto[T"8"lin the electronic version of the paper. 



Target 


SpT 


T c ff 


l°g£ 


ft(0) 




Ae (cr) 


A6 (£) 


Ae o (T eff +500) 


6 O (T eff +500, T efr 


HD 36512 


B0 V 


33700 


4.2 


4.4+1.5 


8.71 


0.10 


0.05 


±0.06 


(8.76, 


8.65) 


HD 37020 


B0.5V 


30500 


4.2 


6.4+1.6 


8.70 


0.10 


0.07 


±0.05 


(8.74, 


8.65) 


HD 36960 


B0.5V 


28900 


3.9 


5.9+0.8 


8.71 


0.10 


0.04 


±0.03 


(8.74, 


8.68) 


HD 37042 


B0.7V 


29700 


4.2 


4.9+1.1 


8.75 


0.08 


0.06 


±0.02 


(8.78, 


8.74) 


HD 36591 


Bl V 


27200 


4.1 


4.5+0.3 


8.71 


0.10 


0.02 


+0.02 


(8.71, 


8.74) 


HD 36959 


Bl V 


25800 


4.2 


2.1+0.4 


8.70 


0.06 


0.02 


+0.05 


(8.66, 


8.76) 


HD 37744 


B1.5 V 


23800 


4.1 


3.6+1.4 


8.70 


0.07 


0.06 


+0.09 


(8.61, 


8.79) 


HD 35299 


B1.5 V 


23700 


4.2 


2.8+0.6 


8.72 


0.07 


0.03 


+0.09 


(8.64, 


8.82) 


HD 36285 


B2V 


20600 


4.0 


5.5+1.5 


8.80 


0.10 


0.06 


+0.13 


(8.67, 


8.92) 


HD 35039 


B2V 


19800 


3.7 


5.3+1.5 


8.79 


0.07 


0.07 


+0.15 


(8.65, 


8.94) 


HD 36629 


B2V 


20000 


4.1 


6.0+1.7 


8.76 


0.10 


0.06 


+0.14 


(8.62, 


8.89) 


HD 36430 


B2V 


18600 


4.1 


6.3+2.2 


8.76 


0.07 


0.08 


+0.13 


(8.64, 


8.90) 


HD 35912 


B2V 


18500 


4.0 


6.3+2.2 


8.79 


0.09 


0.08 


+0.13 


(8.67, 


8.93) 
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Fig. 5 Comparison of our derived stellar parameters, Si and O abundances, with those obtained by ICunha & Lambert) (1992, 
1994 . Results from both studies for each star are connected by a dashed line. Numbers follow the order of stars presented in 
Tables |3] and|4] The size of the uncertainties of the various quantities are shown by the crosses. 



can find slightly different results (but within the errors) when 
analyzing spectra from two different observing campaign^ 



Recently. iPrzybilla. Nieva. & Butler (2008) analyzed a rep- 
resentative sample of unevolved early B-type stars in nearby 
OB associations and the field using a similar technique but 
a different stellar atmosphere code (line blanketed atlas9 lte 
model atmospheres, Kurucz 1993) and nlte line-formation 
calculations (with updated versions of detail and surface, 
Giddings 198 1 ; Butler & Giddings 1985). They have one star in 
common with our work, HD 36591 (HR 1861). We obtain very 
similar results for the stellar parameters, as well as the Si and 
O abundances. They used several T e fj spectroscopic indicators 
(apart from Si ii-iv), finding very good overall agreement. 



Gummersbach et al.l (119981) included a sample of 5 stars 



from Ori OB 1 in their study of the abunda nce gradient of the 



Galaxy. We have three stars in common. iGummersbach et al 



10 The different results may also be caused by actual changes in the 
stellar spectra (e.g. due to binarity). 



used a self-consistent spectroscopic approach. We find similar 
effective temperatures for the two hotter objects (when Siiv 
lines are available), but ~ 0.2 dex lower gravities. For these 
stars (HD 36959 and HD 36960), the derived Si and O abun- 
dances agree (within the errors) with our values. The difference 
in gravity could explain their slightly lower O abundances. 

Interestingly, for the cooler object in common (HD 35039), 
we obtained lower T e fr and logg. It is remarkable that they 
obtained very low O and Si abundances (8.20 and 7.08 dex, 
respectively) for this star. Probably, this result is related to the 
Si ii problem we described in Sect. l3.4l (see also lower panel in 
Fig. [2); they used line Si n4130 to establish T e ff, so obtained 
too high a value (23500 K vs. 19900 K, the value we obtained). 
For this range of stellar parameters, higher T e ff supposes lower 
O abundances, and a change of -3000 K in T e g can perfectly 
explain a 0.5 dex variation in the O abundance. 

Figure H] presents a comparison our r esults with those ob- 
tained by ICunha & Lambert! (119921 1 19941) concerning the stel- 
lar parameters, Si and O abundances. There is a clear discrep- 
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ancy for the majority of the objects, in ter ms not only of ef- 
fective temp eratures, but also of gravities. Cunha & Lambert! 
dl 9921,1 19941) made use of the calibrations of Stromgen photom- 
etry coupled with the fits to the pressure-broadened line wings 
of Hy (from lKuruczlll979|, LTE stellar at mopsheres) to derive 
T e ff and log g. Followi"n glGies & Lambert! dl992l), they adjusted 
the T e ff calibrations by Lester et al.l (119861) and iBalonal (119841) 
upward by 4.2% and 5.2% , respectively. The discrepancy in 
effective temperatures obtained through commonly used pho- 
tometric calibrations and spectroscopic line diag nostic has al- 
ready been pointed out by sev eral authors (e.g. iKilian et all 
19911: iNieva & Przybillal 120081) . Our result once more brings 
out these discrepancies. 

Note also the l arge discrepancy found for the gravity. The 
Cunha & Lambert] values are systemati cally higher (~0.3 dex). 
As shown bv lNieva & Przybillal d2007l) . the use of lte profiles 
for the gravity determination lead to overestimated log g values, 
in particular for hotter stars. In Fig. [3] we show a comparison 
of synthetic hydrogen lines from two fastwind models with the 
observed profiles for HD 36960 (labeled with #3 in Fig. [5] One 
of the models considers the stellar parameters derived from our 
analysis, the other one correspon ds to the T e ff-log g pair pro- 
vided bv lCunha & Lambert! d 19921) . The difference in the wings 
of the lines is clear. 

We also find big discrepancies in the derived Si and O abun- 
dance for most of the stars, with no systematic trend (i.e. dif- 
ferences are found in both positive and negative directions). 
There are several factors to take into account to explain this 
disagreement. One of them is the imprecise determination of 
the stellar parameters from photometric indices. But maybe the 
most important one is the following. Once the stellar param- 
eters were calculated as described above, ICunha & Lambert 
use those values to perform the abundance analysis, comput- 
ing lin e-blanketed lte model atmospheres with atlas6 (IKuruczi 
19791) f or the lte cas e, or using a grid of EWs based on Gold's 
models dGoldll 19841) for the nlte abundance calculations. The 
risk in this strategy is that the stellar abundance analysis is de- 
coupled of the stellar parameter determination, and this could 
produce inconsistencies in the analysis process. Briefly, the cal- 
culation of EWs of those lines used for the abundance determi- 
nation is based on the stellar atmosphere structure defined by 
a stellar atmosphere model computed for a given set of stel- 
lar parameters; on the other hand, the stellar parameters de- 
rived for a given star will depend on the characteristics of the 
stellar atmosphere model we have use Therefore, it may 
be dangerous to do an abundance analysis with a given stellar 
atmosphere code using the stellar parameters obtained from a 
different code or a photometric calibration. This argument is 
crucial in the case of stars with T e ff > 30000 K, where photo- 
metric methods become completely unreliable discriminators 
of temperatures and gravities, because of the insensitivity of 



11 As an example, we want to mention the consequences that in- 
cluding the line blanketing and wind blanketing effects in the stel- 
lar atmosphere models of O and early B-ty pe stars had on the 
SpT-T eff calibrations of these stars (see e.g. Martins et al., 2005; 
iRepolust. Puis, & Herrercll2004h . 



the Rayleigh-Jeans tail of the spectral energy distribution on 
temperature. 

A self-consistent spectroscopic approach allows minimiz- 
ing this problem. In this case, the stellar parameters are deter- 
mined by fitting certain spectroscopic diagnostics, and then the 
same models are used for the abundance analysis. This way 
we are certain that the stellar atmosphere structure used for 
the computation of the abundance diagnostics is coherent with 
the derived stellar parameters for the studied star. To a first or- 
der, when the whole analysis is performed with stellar atmo- 
sphere codes with different characteristics, the derived abun- 
dances should be quite similar (although the derived stellar pa- 
rameters could be somewhat different). This approach can be 
strengthened if multiple independent spectroscopic indicators 
are considered (i.e . Si m- iv, He i-n, C ii-iv; see e.g. this study or 
INieva & Przybillal l2008h . 

5. Chemical composition of B-type stars in the 
OriOBI association 

Figures [6] and [7] again show the derived oxygen and silicon 
abundances in our sample of B-ty pe stars in the OriOBI as - 
sociation, along with the results by ICunha & Lambert! d 19941) . 
This time, the star s are o rdered following t he subgroups sug- 
gested by iBlaauwl dl964h . As discussed in Cunha & Lambert! 
(1994), they found that the stars in the youngest group (Id) 
and some of the stars in subgroup Ic were enriched in oxy- 
gen by about 40% relative to the stars belonging to the older 
subgroups. They interpreted this result as possible proof of 
enrichment of the new generation of stars in the association 
with the products from supernovae ejecta from the older sub- 
groups. In addition, they found features of this enrichment in 
silicon (also expected from the triggered star formation sce- 
nario). In contrast, we did not find any systematic difference 
between the O and Si abundances in stars from the various as- 
sociations. In fact, our results indicate that the B-type stars in 
the Ori OB 1 association are chemically homogeneous (at least 
in terms of oxygen and silicon), having a dispersion in abun- 
dances (0.04 and 0.03 dex, respectively) smaller than the in- 
trinsic uncertainties of the derived abundances (0.10 and 0.08 
dex, respectively). The mean abundances are e(0) = 8.73 dex 
and e(Si) = 7.5 1 dex. These values agree with those obtained by 
IPrzybilla. Nieva. & Butlerl (12008b for their sample of six stars 
in the solar neighborhood. 



5.1. Comparison with the Sun 

We also include in Fig. [6] the resulting solar oxygen abun- 
dances appearing in the literature in the past 4 years (during 
the so-called solar oxygen crisi^H epoch, not yet finished). 
Various authors have published values for the solar oxygen 
abundance based on improved model atmospheres (either ID 
and 3D), line formation codes, atomic and molecular data, and 
detai led treatment of blends. I n the p l ot, we present results 
by Centeno & Socas-Navarrol ( 2008I) . Melendez & Asplund 
(2008), Avresl d2008l) . ICaffau et all d2008l) . lAllende Prieto 



12 lAvres et al. (2006) 
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Fig. 6 Oxygen abund ances derived f or our sample of early B-type stars in the Orion association, and comparison with the previous 
results from lCunha & Lambert! i 19941) . Vertical lines separates stars from the various subgroups. Solid and dashed horizontal lines 
represent the mean value and the dispersion (lcr) of our results; dot-dashed horizontal lines indicate the characteristic intrinsic 
uncertainty of the derived abundances (accounting for unc ertainties in stellar p arameters, microturbulence, and the line-by-line 
abundance dispersion). The oxygen abundance derived by Esteban et al.l ( 120041) for M42, along with the values determined for 
the Solar abundance in the past 4 years (the epoch of the "solar crisis"), are also presented for comparison. 
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Fig. 7 As Fig.[2]but for silicon. This time, the referen ce for the nebular abundance is Rubin et al. d 1993b and Garnett et aL , 
(first and second column, see text). Solar values from Asplund. Grevesse. & Sauval ( 20051) and Grevesse & SauvaT i 1998 ) 



(1995) 



(20 081 ISocas-Navarro & Nortonl d2007l). lAvres et al] (I2006I) . 
lAsplund et al.1 d2004l) . and lGrevesse & Sauvall(ll998b . This last 
value (marked with an internal black dot) was considered as the 
standard solar O abundance until a f ew years ago. The derived 
solar values range from 8.63 dex ( Socas-Navarro & Norton , 
2007) to 8.86 dex dCenteno & Socas-Navarrol 120081) . TheO 



of all these values. In view of the present-day results, the only 
thing we can say is that oxygen abundances in the Sun and 
B-type stars in the solar vicinity are the same within the un- 
certainties. However, we consider it too premature to draw any 
firm conclusion or hypothesis about the chemical evolution of 
the local interstellar medium during the lifetime of the Sun. 



abundances in our sample of stars in Ori OBI lie in the middle 
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Figure [7] also shows a comparison of our derived Si abun- 
dances witlrtlie^oja^ralue.^^only present here the old abun- 
dance (Gr evesse & Sa uvaill998l) and the new value proposed 
by lAsplund. Grevesse. & Sauval d2005l) . In contrast to previ- 
ous results, the mean Si abundance in B-type stars in Ori OB 1 
is very close to t he Solar value (similar to other B- type stars in 
the solar vicinity IPrzybilla. Nieva. & Bufleru2008l) . 



5.2. Comparison with the Orion nebula 

The most recent and detailed anal ysis of the opti c al spe ctrum 
of the Orion nebula was done by Esteban et al.1 (120041) . The 
high-quality UVES@VLT spectrum they used allowed them to 
derive the oxygen gas phase abundance of the nebula by us- 
ing collisionally excited lines (CELs) and recombination lines 
(RLs). The final value they propose is 8.65+0.03 dex. This 
value was calculated assuming the ionic abundances given by 
2 + RLs and + CEL (plus a t 2 =O.O220 Note, however, 
the abundance given by the 2 + and + CELs (and t 2 =0) 
is 8.51+0.03 dex. These values have been included in Fig. [6] 
for comparison with the early B-type stellar abundances pre- 
sented in this study. The mean value of the derived stellar oxy- 
gen abundances is 0.25 and 0.11 dex higher than the nebular 
abundances given by the CELs and RLs, respectively. 

Although the mean stellar abundance seems to agree better 
with the one given by the faint recombination lines, we have to 
consider that the analysis of the nebular emission line spectrum 
can only provide abundances for the ionized gas phase of the 
ISM. This means a lower limit to the actual ISM abundance, 
since part of the oxygen can be depleted, forming part of the 
dust grains. 

Silicon is one of the elements expected to be more dep leted 
onto dust grains, along with Mg and Fe dDrainel 120031) . We 
can therefore compare our derived Si abundances with those 
obtained from the study of the emission line spectrum of the 
Orion nebula to try to find some clues about the amount of 
oxygen depleted. One should note that the determination of Si 
nebular abundances is not straightforward and must be based 
on results from photoionization models and certain ionization 
correction factors (ICFs). The only determinations of Si abun- 
dance in the Orion nebula we found in the literature are those 



by iRubinetal J (119931) and iGarnett et all (119951) . Both studies 
used the same observations of the FUV Si ra] /LU883, 1892, 
and C m] /L11907, 1909 lines to estimate the Si abundance, and 
obtained 6.65 and 6.58 dex, respec tively. They comp uted pho- 
toionization models to obtain Si/H (IRubin et all 1 19931) and the 
ICF needed to transform Si 2+ / C 2+ into Si/C henc e the C/O 



and O/H ratios to derive Si/H (IGarnett et all 119951) . In addi 



Garnett et alj ( 19951) discus s how the effect of assuming 



tion, 

a t 2 =0.04 dPeimbertet all 1 1993b in the calculation of the C/O 



and O/H ratios could affect the derived Si abundance. They ob- 
tained log (Si/H)+12 = 7.14 dex for this case. These three val- 



13 The t 2 parameter was introduced by IPeimberj i 19671) to account 
for temperature fluctuations in ionized nebulae. The temperature fluc- 
tuation scenario has be en proposed to explain the CEL s vs. RL abun- 
dance discrepancy (see lGarcia-Roias & Estebanl 2007, and references 
therein). 



ues of the Si abundance are included in Fig [7] for comparison 
with the stellar abundances. The difference between stellar and 
gas phase Si abundances are ~1 dex, or ~0.3 dex, depending 
on the nebular abundances we trust more. 

Since a detailed comparison of stellar and nebular abun- 
dances in Ori OBI and the Orion nebula within a dust-depletion 
scenario requires more extended study, which is beyond the 
scope of this paper, we decide to present a more detailed dis- 
cussion in a separate paper (in preparation). 

6. Summary and conclusion 

In this work, we applied a self-consistent spectroscopic ap- 
proach to determining of the stellar parameters, the Si and O 
abundances of a sample of early-B type stars from the various 
subgroups of the Ori OB 1 association. We made use of a high- 
quality spectroscopic data set, obtained with fies@not, and the 
modern nlte, line-blanketed, spherically extended stellar atmo- 
sphere COde FASTWIND. 

We developed several IDL programs to automatically iden- 
tify and measure the EW of metal lines in high-resolution spec- 
tra. We also constructed a grid of fastwind HHeSiO models op- 
timized for the analysis of early B-type, main sequence stars. 

The availability of a large number of Si ii-m-iv lines in the 
fies@not spectra allowed us to obtain the stellar parameters 
with high accuracy and detect some problems related to some 
Si lines commonly used for the stellar parameter and Si abun- 
dance determination. Once these problems were accounted for, 
a high degree of homogeneity was found for Si abundances in 
the analyzed sample of stars. 

The oxygen abundance analysis also result in a small 
dispersion of abundances, in contrast to previous determina- 
tions. The mean oxygen and silicon abundances agree with 
those resulting from a similar analysis of a representative 
sample of unevolved early B-typ e stars in the solar vicinity 
dPrzvbilla. Nieva. & Butlen, (2008). Both results indicate that 
abundances derived from these stellar objects are more homo- 
geneous and metal-rich than previously though. 

We also compared the O and Si stellar abundance in 
Ori OB 1 with those obtained for the Sun during the epoch of 
the "solar crisis". The O abundances in our sample of stars in 
Ori OBI lie in the middle of all these values. In view of the 
present-day results, the only thing we can say is that oxygen 
abundances in the Sun and B-type stars in the solar vicinity 
are the same within the uncertainties. However, we consider it 
too premature to draw any firm conclusion or hypothesis about 
the chemical evolution of the local interstellar medium during 
the lifetime of the Sun. Silicon abundances are also very sim- 
ilar (contrary to what was previously found from the study of 
B-type stars). 

Finally, we compared the stellar abundances with those de- 
rived from the study of the emission line spectrum of the Orion 
nebula. In a forthcoming paper we will present a more detailed 
discussion accounting for the possible depletion of O, Si, Mg, 
and Fe into dust. 

This work points out one more time the importance of self- 
consistent spectroscopic abundance analyses for determining 
of the chemical composition of the photospheres of OB-type 
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stars. Photometric T e ff diagnostics must be treated with cau- 
tion in the context of abundance analyses of these type of ob- 
jects. It is also dangerous to combine the stellar parameters 
determined with any given code and the abundance analysis 
performed with a different code. Systematic errors inherent in 
those techniques, or possible biases between the various codes, 
can lead to incorrect abundances. 
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Table 5 List of log gf values for the Si ii-iv lines considered in 
this stud y (from Atomic Line Li st v2.05). For the On lines we 
refer to Simon- Diaz etaD d2006l) 



Ion 


A A 


loe ef 

o 6./ 


Configuration 


Term 


ji 11 


JoJO.UZ 




1c In 2 1c 2 t\ SL\ A.W 


in ?p n 

_ \J- —V\) 


Si ii 


3862.60 


-0.7130 


3s 3p 2 -3s 2 (IS) 4p 


2D-2Po 


Sin 


4128.10 


0.3646 


3s 2 (IS) 3d-3s 2 (IS) 4f 


2D-2Fo 


Sin 


4130.90 


0.5192 


3s 2 (IS) 3d-3s 2 (IS) 4f 


2D-2Fo 


Si ii 


5041.02 


0.2543 


3s 2 (IS) 4p-3s 2 (IS) 4d 


2Po-2D 


Sin 


5055.98 


0.5083 


3s 2 (lS)4p-3s 2 (lS)4d 


2Po-2D 


Sin 


5056.32 


-0.4460 


3s 2 (lS)4p-3s 2 (lS)4d 


2Po-2D 


Sin 


6347.11 


0.1819 


3s 2 (IS) 4s-3s 2 (lS)4p 


2S-2Po 


Sin 


6371.37 


-0.1208 


3s 2 (IS) 4s-3s 2 (lS)4p 


2S-2Po 


Sim 


4552.62 


0.2828 


3s 4s-3s 4p 


3S-3Po 


Sim 


4567.84 


0.0595 


3s 4s-3s 4p 


3S-3Po 


Sim 


4574.76 


-0.4183 


3s 4s-3s 4p 


3S-3Po 


Sim 


5739.73 


-0.1032 


3s 4s-3s 4p 


lS-lPo 


Si iv 


4088.86 


0.1984 


2p 6 4s-2p 6 4p 


2S-2Po 


Si iv 


4116.10 


-0.1055 


2p 6 4s-2p 6 4p 


2S-2Po 


Si iv 


4212.40 


0.3840 


2p 6 5d-2p 6 6f 


2D-2Fo 


Si iv 


4212.41 


0.5601 


2p 6 5d-2p 6 6f 


2D-2Fo 


Si iv 


4631.24 


1.2158 


2p 6 5f-2p 6 6g 


2Fo-2G 



S. Simon-Diaz: O and Si abundances in B-type stars in Ori OB 1, Online Material p 3 



Table 6 Results from the abundance analysis of HD 365 12 (BO V) 
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Table 7 Results from the abundance analysis of HD 37020 (B0.5 V) 
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10 


8.73 


0.15 


OH4891 


40 


8 


8.84 


0.19 


OII4906 


48 


9 


8.63 


0.18 


OII4941 


48 


6 


8.57 


0.12 


OII4943 


65 


8 


8.53 


0.13 








= 8.70 


Ae o (o-) = 0.10 




Af,(0) = 1.6 


=> 


Af O (f,) = 0.07 



eo=8-71 Ae o (o-) = 0.10 
Af,(0) = 1.5 => Ae o (f.) = 0.05 
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Table 8 Results from the abundance analysis of HD 36960 (B0.5 V) 



HD 36960 


T c rr = 


28900 K, log 


g = 3.9 dex 


f,(Si)=5.4 


Line 


EW 


AEW 


^Si 






(mA) 


(mA) 




(dex) 


SiIII4552 


211 


5 


7.53 


0.06 


SiIII4567 


179 


5 


7.51 


0.06 


SiIII4574 


116 


5 


7.52 


0.07 


SiIV4089 


160 


11 


7.55 


0.15 


SiIV4116 


127 


9 


7.52 


0.14 


SiIV4212 


28 


4 


7.50 


0.14 


SiIV4631 


37 


3 


7.54 


0.08 








£ S i = 7-53 


Ae S i(<x) -0.02 




Af,(Si) = 0.6 


— \ 


A& ■({.") -0 06 




T e ff - 


28900 K, log 


g = 3.9 dex 


ft(P)=5.9 


OII3945 


68 


6 


8.78 


0.11 


OII3954 


96 


7 


8.71 


0.10 


OII3982 


61 


6 


8.63 


0.11 


OII4317 


125 


7 


8.71 


0.08 


OII4319 


125 


6 


8.78 


0.08 


OII4366 


113 


5 


8.66 


0.07 


OII4414 


165 


5 


8.66 


0.06 


KJLl^H I 


141 


5 


8.75 


U.UO 


OII4452 


50 


5 


8.54 


0.10 


OII4661 


127 


5 


8.87 


0.08 


OII4673 


46 


5 


8.85 


0.12 


OII4676 


101 


5 


8.66 


0.08 


OII4696 


30 


3 


8.78 


0.09 


OII6641 


47 


6 


8.81 


0.12 


OII6721 


76 


9 


8.77 


0.12 


OII4078 


60 


4 


8.81 


0.09 


OII4891 


34 


10 


8.64 


0.27 


OII4906 


61 


6 


8.71 


0.11 


OII4941 


54 


10 


8.54 


0.19 


OII4943 


75 


4 


8.54 


0.06 


OII4956 


23 


4 


8.72 


0.14 






eo = 8-71 


Ae o (cr) = 0.10 




Afi(O) = 0.8 


=> 


Ae o (?,) = 0.04 



Table 9 Results from the abundance analysis of HD 37042 (B0.7 V) 



HD 37042 


T c( r = 


29700 K, log g 


= 4.2 dex 


ft(Si)=1.4 


Line 


EW 


AEW 


esi 






(mA) 


(mA) 


(dex) 


(dex) 


SiIII4552 


149 


2 


7.54 


0.04 


SiIII4567 


133 


6 


7.59 


0.11 


SiIII4574 


88 


4 


7.56 


0.08 


SiIII5739 


99 


7 


7.50 


0.13 


SiIV4116 


95 


10 


7.53 


0.21 


SUV4212 


24 


5 


7.54 


0.22 


SiIV4631 


35 


10 


7.55 


0.30 






e si = 7.55 


Aesi(o-) = 0.03 




Af t (Si) = 0.3 


=> 


Aesi(f,) = 0.04 


HD 37042 


Tcir = 


29700 K, log g 


= 4.2 dex 


ft(0)=4.9 


OII3945 


68 


6 


8.87 


0.11 


OII3954 


93 


7 


8.79 


0.12 


OII4317 


116 


10 


8.75 


0.13 


OII4319 


116 


11 


8.84 


0.16 


OII4366 


108 


9 


8.76 


0.14 


OII4414 


143 


10 


8.64 


0.13 


OII4416 


127 


7 


8.78 


0.10 


OII4452 


64 


6 


8.82 


0.12 


OII4641 


142 


7 


8.81 


0.11 


OII4661 


108 


5 


8.79 


0.09 


OII4676 


95 


6 


8.71 


0.11 


OII4696 


30 


5 


8.82 


0.16 


OII4072 


110 


10 


8.62 


0.18 


OII4076 


128 


10 


8.63 


0.17 


OII4078 


57 


8 


8.84 


0.19 


OII4891 


39 


5 


8.80 


0.13 


OII4906 


57 


5 


8.75 


0.10 


OII4941 


53 


4 


8.63 


0.08 


OII4943 


73 


4 


8.63 


0.07 


OII4956 


23 


8 


8.79 


0.29 








= 8.75 


Ae o (o-) = 0.08 




Af,(0) = l.l 


=> 


Ae o (f,) = 0.06 
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Table 1 Results from the abundance analysis of HD 36591 (Bl V) 



HD 36591 


Tefl = 


27200 K, log 


g = 4. 1 dex 


ft(Si)=1.3 


Line 


EW 


AEW 








(mA) 


(mA) 


(dex) 


(dex) 


SiIII4552 


158 


2 


7.52 


0.03 


SiIII4567 


145 


2 


7.62 


0.04 


SiIII4574 


104 


2 


7.65 


0.04 


SiIII5739 


114 


2 


7.56 


0.03 


SiIV4089 


80 


2 


7.54 


0.05 


SiIV4116 


61 


2 


7.48 


0.06 








e si = 7.53 


Ae si (cr) = 0.06 




Af,(Si) = 0.3 


=> 


Ae si (f t ) = 0.03 


HD 36591 


T c ,r = 


27200 K, log 


g = 4. 1 dex 


f,(0)=4.5 


OII3945 


69 


2 


8.80 


0.04 


OII3954 


86 


2 


8.64 


0.04 


OII3982 


62 


2 


8.64 


0.04 


OII4317 


125 


5 


8.77 


0.06 


OII4319 


114 


2 


8.76 


0.03 


OII4366 


102 


4 


8.62 


0.07 


OII4414 


144 


2 


8.61 


0.03 


OII4416 


124 


8 


8.70 


0.12 


OII4452 


53 


3 


8.59 


0.07 


OII4638 


109 


2 


8.94 


0.04 


OII4641 


145 


2 


8.84 


0.03 


OII4650 


159 


4 


8.62 


0.06 


OII4661 


113 


2 


8.83 


0.03 


OH4673 


44 


5 


8.81 


0.13 


OH4676 


97 


4 


8.70 


0.07 


OH4696 


30 


3 


8.74 


0.10 


OII6641 


48 


2 


8.90 


0.04 


OII6721 


69 


2 


8.79 


0.03 


OII4072 


113 


7 


8.65 


0.13 


OII4076 


133 


4 


8.70 


0.07 


OII4078 


56 


3 


8.77 


0.08 


OII4086 


68 


6 


8.83 


0.14 


OII4891 


34 


3 


8.67 


0.09 


OII4906 


53 


3 


8.67 


0.07 


OII4941 


48 


2 


8.54 


0.04 


OII4943 


68 


3 


8.57 


0.06 


OII4956 


24 


3 


8.77 


0.11 



Table 11 Results from the abundance analysis of HD 36959 (B1V) 



HD 36959 


T c ,r = 


25800K.log 


g = 4.2 dex 


f t (Si)=0.0 


Line 


EW 


AEW 








(mA) 


(mA) 


(dex) 


(dex) 


SiII4128 


16 


3 


7.52 


0.12 


SiI14130 


17 


3 


7.44 


0.12 


SiIII4552 


142 


1 


7.48 


0.02 


SiIII4567 


124 


1 


7.48 


0.02 


SiIII4574 


90 


1 


7.54 


0.02 


SiIII5739 


99 


2 


7.51 


0.04 


SiIV4089 


45 


1 


7.45 


0.04 


SiIV41 16 


36 


3 


7.47 


0.12 








e si = 7.49 


Ae si (cr) = 0.03 




Aft(Si) = 0.5 


=> 


Ae si (£) = 0.07 


HD 36959 


T cf r = 


25800 K, log 


g = 4.2 dex 


f,(0)=2.1 


OII3945 


55 


2 


8.78 


0.05 


OII3954 


73 


2 


8.69 


0.04 


OII3982 


50 


2 


8.63 


0.05 


OII4317 


98 


6 


8.76 


0.10 


OII4319 


86 


3 


8.69 


0.06 


OII4366 


81 


5 


8.62 


0.10 


OII4414 


116 


2 


8.64 


0.03 


OII4416 


99 


3 


8.72 


0.05 


OII4452 


45 


2 


8.62 


0.05 


OII4641 


113 


2 


8.75 


0.04 


OII4650 


133 


5 


8.65 


0.08 


OII4661 


89 


2 


8.76 


0.04 


OII4673 


35 


2 


8.72 


0.07 


OII4676 


75 


3 


8.61 


0.07 


OII4696 


25 


2 


8.69 


0.08 


OII6641 


30 


2 


8.78 


0.06 


OII6721 


48 


2 


8.75 


0.04 


OII4072 


94 


6 


8.72 


0.13 


OII4076 


108 


6 


8.73 


0.12 


OII4078 


46 


3 


8.80 


0.09 


OII4086 


48 


3 


8.70 


0.09 


OII4891 


26 


2 


8.67 


0.07 


OII4906 


40 


2 


8.67 


0.06 


OII4941 


39 


2 


8.62 


0.06 


OII4943 


55 


2 


8.66 


0.05 



eo=8-71 Ae o (o-) = 0.10 
A£(O) = 0.3 Ae o (£) = 0.04 



e o =8-70 Ae o (o-) = 0.02 
Af,(O) = 0.4 =• Ae o (f,) = 0.06 
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Table 1 2 Results from the abundance analysis of HD 37744 (B 1 .5 V) 



HD 37744 


Toff = 


23800 K, log 


g = 4. 1 dex 


ft(Si)=0.5 


Line 


EW 


AEW 








(mA) 


(mA) 




taexj 


SiII3856 


20 


4 


7.49 


0.22 


SiII6371 


18 


6 


7.57 


0.22 


SiIII4552 


147 


3 


7.57 


0.05 


SiIII4567 


120 


4 


7.49 


0.07 


SiIII4574 


88 


4 


7.59 


0.09 


SiIII5739 


86 


4 


7.45 


0.08 


SiIV4116 


23 


4 


7.56 


0.26 








£ S i = 7-53 


Ae S i(<x) -0.05 




Af t (Si) = 0.5 


— \ 


A& ■({.") -0 04 


un inn a A 


Teff = 


23800 K, log 


g = 4.1 dex 


ft(O)=3.0 


OII3945 


45 


5 


8.73 


0.13 


OII3954 


66 


7 


8.72 


0.15 


OH4317 


72 


7 


8.60 


0.13 


OH4319 


73 


8 


8.67 


0.16 


OII4366 


67 


6 


8.56 


0.13 


OH4414 


100 


14 


8.62 


0.23 


OH4416 


76 


7 


8.59 


0.14 




32 


6 


8.51 


n on 
u.zu 


OII4641 


99 


5 


8.76 


0.09 


OII4650 


123 


8 


8.71 


0.13 


OII4661 


79 


5 


8.81 


0.10 


OII4673 


25 


5 


8.67 


0.20 


OH4676 


71 


6 


8.76 


0.13 


OII4696 


20 


5 


8.73 


0.23 


OII6721 


33 


7 


8.74 


0.19 


OII4069 


136 


13 


8.83 


0.18 


OII4078 


36 


6 


8.71 


0.21 


OII4906 


32 


4 


8.69 


0.13 


OH4941 


30 


5 


8.62 


0.17 


OII4943 


45 


5 


8.69 


0.13 


OII4956 


14 


3 


8.79 


0.17 






e = 8-70 


Ae (cr) = 0.07 




Af,(0) = 1.4 


=> 


A £o (?,) = 0.06 



Table 13 Results from the abundance analysis of HD 35299 (B 1.5 V) 



HD 35299 


T cf r = 


23700 K, log g = 4.2 dex 


ft(Si)=0.5 


Line 


EW 


AEW 


<5i 






(mA) 


(mA) 


(dex) 


(dex) 


SUI3856 


27 


5 


7.65 


0.24 


SiII3862 


15 


2 


7.44 


0.13 


SiII6371 


23 


4 


7.64 


0.13 


SiIII4552 


133 


5 


7.45 


0.08 


SiIII4567 


110 


2 


7.40 


0.04 


SiIII4574 


79 


2 


7.50 


0.05 


SiIII5739 


81 


2 


7.47 


0.04 


SiIV4089 


26 


1 


7.59 


0.06 


SiIV4116 


18 


2 


7.48 


0.15 








e si = 7.50 


Ae si (o-) = 0.08 




Af,(Si) = 0.5 


=> 


Ae si (f,) = 0.02 


HD 35299 


T c ,r = 


23700 K, log g = 4.2 dex 


ft(0)=2.8 


OII3945 


44 


2 


8.81 


0.06 


OII3954 


59 


2 


8.72 


0.05 


OII3982 


38 


2 


8.62 


0.06 


OII4317 


78 


5 


8.80 


0.09 


OII4319 


66 


3 


8.69 


0.07 


OII4366 


61 


3 


8.59 


0.07 


OII4414 


95 


5 


8.70 


0.09 


OII4416 


80 


2 


8.77 


0.04 


OII4452 


32 


2 


8.61 


0.07 


OII4641 


91 


5 


8.79 


0.10 


OII4650 


112 


5 


8.73 


0.09 


OII4661 


68 


3 


8.77 


0.07 


OII4673 


25 


2 


8.76 


0.08 


OII4676 


59 


2 


8.68 


0.05 


OII4696 


17 


2 


8.70 


0.10 


OII6641 


18 


2 


8.82 


0.09 


OII6721 


29 


3 


8.77 


0.09 


OII4072 


82 


3 


8.81 


0.07 


OII4076 


94 


4 


8.79 


0.08 


OII4078 


35 


2 


8.79 


0.07 


OII4086 


36 


2 


8.67 


0.07 


OII4891 


16 


3 


8.64 


0.16 


OII4906 


28 


2 


8.70 


0.07 


OII4941 


25 


2 


8.59 


0.08 


OII4943 


40 


2 


8.71 


0.06 


OII4956 


12 


2 


8.78 


0.13 



eo=8-72 Ae o (o-) = 0.07 
Af,(O) = 0.6 Ae o tf,) = 0.03 
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Table 1 4 Results from the abundance analysis of HD 36285 (B2 V) 



HD 36285 


Toff = 


20600 K, log 


g = 4.0 dex 


ft(Si)=1.7 


Line 


EW 


AEW 


^Si 


Aesi 




(mA) 


(mA) 


(dex) 


(dex) 


SiII3856 


41 


2 


7.38 


0.07 


SiII3862 


37 


5 


7.52 


0.19 


SiII6347 


61 


4 


7.60 


0.09 


SiII637 1 


45 


3 


7.53 


0.08 


SiIII4552 


110 


1 


7.49 


0.02 


SiIII4567 


88 


2 


7.44 


0.04 


SiIII4574 


59 


1 


7.51 


0.03 


SiIII5739 


53 


2 


7.46 


0.05 








e Si = 7-49 


Ae si (o-) = 0.06 




Af t (Si) = 0.5 


=> 


Ae si (f t ) = 0.05 


HD 36285 


Teff = 


20600 K, log 


g = 4.0 dex 


<f,(0)=5.5 


0113945 


26 


3 


8.79 


0.12 


0113954 


41 


3 


8.79 


0.09 


0114317 


46 


5 


8.74 


0.12 


0114319 


42 


3 


8.69 


0.08 


0114366 


42 


4 


8.65 


0.11 


0114414 


63 


3 


8.70 


0.06 


/""ITT A A I C 
0114416 


50 


3 


8.74 


0.07 


OII4452 


26 


3 


8.88 


0.12 


OII4638 


46 


3 


8.94 


0.08 


OII4650 


78 


4 


8.70 


0.07 


OII4661 


46 


2 


8.80 


0.05 


OII4673 


12 


3 


8.68 


0.21 


OH4676 


36 


3 


8.66 


0.09 


OH6721 


15 


4 


8.99 


0.23 


OII4069 


81 


6 


8.86 


0.10 


OII4072 


55 


4 


8.80 


0.11 


OII4078 


24 


4 


8.91 


0.19 


OII4086 


24 


10 


8.78 


0.48 


OII4906 


16 


3 


8.85 


0.17 


OII4941 


16 


5 


8.84 


0.29 


OII4943 


25 


3 


8.91 


0.12 






e = 8-80 


Ae o (<x) = 0.10 




Af,(0) = 1.5 


=> 


Ae o (ft) = 0.06 



Table 15 Results from the abundance analysis of HD 35039 (B2 V) 



HD 35039 


T c ,r = 


19800 K,logg = 3.7 dex 


ft(Si)=3.3 


Line 


EW 


AEW 


<3i 


Aesi 




(mA) 


(mA) 


(dex) 


(dex) 


SiII3856 


61 


2 


7.63 


0.06 


SiII3862 


45 


5 


7.54 


0.15 


SiII6347 


75 


10 


7.51 


0.17 


SiII6371 


55 


2 


7.46 


0.04 


SiIII4552 


125 


10 


7.49 


0.15 


SiIII4567 


103 


5 


7.47 


0.09 


SiIII4574 


67 


5 


7.52 


0.11 


SiIII5739 


69 


5 


7.58 


0.10 








<=si = 7.52 


Ae si (o-) = 0.06 




Af,(Si) = 1.0 


=> 


Ae si (£) = 0.08 


HD 35039 


T cf r = 


19800 K, logg = 3.7 dex 


ft(0)=5.3 


OII3945 


31 


2 


8.91 


0.07 


OII3954 


41 


2 


8.79 


0.06 


OII3982 


25 


2 


8.69 


0.08 


01143 1 / 


50 


6 


8.80 


0.14 


OII4319 


45 


4 


8.74 


0.11 


OH4414 


64 


2 


8.72 


0.04 


OH4416 


54 


4 


8.82 


0.10 


OII4641 


64 


3 


8.80 


0.07 


OII4650 


82 


2 


8.74 


0.04 


OII4661 


46 


3 


8.79 


0.08 


OII4673 


15 


4 


8.79 


0.25 


OH4676 


39 


2 


8.71 


0.06 


OII4069 


80 


3 


8.88 


0.05 


OII4072 


59 


3 


8.91 


0.08 


OII4078 


21 


4 


8.84 


0.21 


OII4086 


25 


4 


8.84 


0.18 


OII4941 


12 


3 


8.68 


0.21 


OII4943 


21 


3 


8.81 


0.14 








eo = 8.79 


Ae o (o-) = 0.07 




Af,(0) = 1.5 


=> 


Ae o (f,) = 0.07 
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Table 1 6 Results from the abundance analysis of HD 36629 (B2 V) 



HD 36629 


T c rr = 


20000 K, log 


g = 4. 1 dex 


ft(Si)=1.0 


Line 


EW 


AEW 


^Si 


Aesi 




(mA) 


(mA) 


(dex) 


(dex) 


SiII3856 


52 


2 


7.58 


0.07 


SiII3862 


42 


2 


7.55 


0.08 


SiII6347 


65 


3 


7.62 


0.07 


SiII637 1 


48 


2 


7.53 


0.06 


SiIII4552 


90 


1 


7.52 


0.02 


SiIII4567 


73 


1 


7.49 


0.02 


SiIII4574 


46 


1 


7.51 


0.03 


SiIII5739 


42 


2 


7.53 


0.06 








e S i = 7.54 


Ae si (o-) = 0.04 




Aft(Si) = 0.5 


=> 


Ae si (f t ) = 0.05 


HD 36629 


Tefl = 


20000 K, log 


g = 4. 1 dex 


£,(O)=6.0 


OII3945 


24 


2 


8.91 


0.09 


OII3954 


28 


2 


8.68 


0.07 


OII4317 


37 


3 


8.77 


0.09 


OII4319 


30 


2 


8.63 


0.07 


0114366 


30 


2 


8.59 


0.07 


OII4414 


50 


2 


8.69 


0.05 


OII4416 


40 


2 


8.75 


0.06 


OII4452 


16 


1 


8.74 


0.06 


OII4638 


36 


1 


8.95 


0.03 


OII4641 


41 


1 


8.63 


0.03 


OII4650 


59 


3 


8.65 


0.07 


OII4661 


32 


2 


8.72 


0.07 


OII4673 


11 


3 


8.82 


0.23 


OII4676 


27 


2 


8.65 


0.07 


OII4069 


64 


3 


8.88 


0.06 


OII4072 


45 


2 


8.83 


0.06 


OII4076 


55 


3 


8.85 


0.08 


OII4078 


15 


2 


8.79 


0.13 


OII4943 


16 


3 


8.83 


0.17 






to = 8.76 


Ae o (o-) = 0.10 




Af,(0) = 1.7 


=> 


A£ O (f t ) = 0.06 



Table 18 Results from the abundance analysis of HD 35912 (B2V) 



HD 35912 


T cf r = 


1 8500 K, log g = 4.0 dex 


f,(Si)=3.2 


Line 


EW 


AEW e S \ 


Atsi 




(mA) 


(mA) (dex) 


(dex) 


SiII3856 


65 


2 7.35 


0.05 


SiII3862 


62 


4 7.53 


0.11 


SH16347 


99 


2 7.60 


0.03 


SM6371 


71 


3 7.49 


0.06 


SiIII4552 


77 


2 7.48 


0.04 


SiIII4567 


60 


2 7.45 


0.05 


SiIII4574 


35 


1 7.47 


0.03 


SiIII5739 


30 


2 7.51 


0.07 






e S i = 7-48 


Ae si (o-) = 0.07 




Af t (Si) = 0.5 => 


A<r si (f,) = 0.04 


HD 35912 


T cf r = 


1 8500 K, log g = 4.0 dex 


«0)=6.3 


OII3945 


15 


5 8.88 


0.32 


OII3954 


22 


4 8.79 


0.18 


OH4317 


25 


5 8.82 


0.20 


OII4319 


21 


3 8.70 


0.14 


OII4366 


22 


5 8.70 


0.23 


OH4414 


34 


4 8.72 


0.13 


OH4416 


29 


3 8.85 


0.11 


OII4650 


42 


2 8.70 


0.06 


OII4661 


25 


3 8.86 


0.13 


OII4676 


17 


4 8.65 


0.22 


OII4069 


44 


3 8.91 


0.08 


OII4072 


33 


3 8.91 


0.12 






£o = 8-79 


Aeo(o-) = 0.09 




Af,(0) = 2.2 => 


Ae o (f,) = 0.08 



Table 1 7 Results from the abundance analysis of HD 36430 (B2 V) 



HD 36430 




18400 K, log 


g = 4. 1 dex 


ft(Si)=3.5 


Line 


EW 


AEW 


ESi 


Ae S i 




(mA) 


(mA) 


(dex) 


(dex) 


SiII3862 


61 


5 


7.43 


0.13 


SiII6347 


100 


10 


7.57 


0.16 


SiII6371 


74 


6 


7.49 


0.12 


SiIII4552 


77 


3 


7.55 


0.06 


SiIII4567 


53 


2 


7.40 


0.05 


SiIII4574 


28 


5 


7.37 


0.19 


SiIII5739 


25 


2 


7.47 


0.08 








e si = 7.47 


Ae si (o-) = 0.08 




Af t (Si) = 1.0 


=> 


Af Si (£) = 0.07 



HD 36430 


Tel 


= 18400 K, logg 


= 4. 1 dex 


f,(0)=6.3 


OII3954 


19 


5 


8.77 


0.26 


OII4317 


20 


5 


8.75 


0.24 


OII4319 


20 


5 


8.76 


0.25 


OII4366 


18 


4 


8.64 


0.21 


OII4414 


33 


3 


8.78 


0.10 


OII4416 


26 


7 


8.85 


0.29 


OII4641 


27 


5 


8.74 


0.20 


OII4650 


36 


2 


8.66 


0.06 


OII4661 


18 


3 


8.72 


0.16 


OII4676 


18 


4 


8.78 


0.21 


OII4069 


37 


3 


8.86 


0.09 


OII4072 


27 


4 


8.82 


0.18 








= 8.76 


Aeo(o-) = 0.07 




Aft(0) = 2.2 


=> 


Aeo(f t ) = 0.08 
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Fig. 8 The complete atlas of fies@not spectra (part 1 of 3). The Si ii-iv and O n lines used for the abundance analysis 
as solid and dashed vertical lines, respectively. 



are indicated 
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Fig. 10 The complete atlas of fies@not spectra (part 3 of 3). The Sin-iv and On lines used for the abundance analysis are 
indicated as solid and dashed vertical lines, respectively. 



